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Abstract 

We find a large class of quantum gauge models with massless fermions where the 
coupling to the gauge fields is not chirally symmetric and which nevertheless do not 
suffer from gauge anomalies. To be specific we study two dimensional Abelian models 
in the Hamiltonian framework which can be constructed and solved by standard tech- 
niques. The general model describes A'p photon fields and Np flavors of Dirac fermions 
with 2N-pNp different coupling constants i.e. the chiral component of each fermion can 
be coupled to the gauge fields differently. We construct these models and find conditions 
so that no gauge anomaly appears. If these conditions hold it is possible to construct 
and solve the model explicitly, so that gauge- and Lorentz invariance are manifest. 



1. Introduction 

In this paper we discuss a large class of gauge theory models where the gauge fields are 
coupled to A'^p Dirac fermions with different coupling constants (= charges) e"^ where the 
index a = +,— distinguishes the chiral components of the fermions, j = 1,2,... ,Np is 
a fermion flavor index, and a = 1,2,... Ap labels different gauge fields. As we shall see, 
the existence of large gauge transformations implies that the charges have to be quantized 
i.e. e"^ = n°^e° (no summation) with n"^ integers. We say that such a model is chirally 
symmetric if for all a and j, _,_ = ''^"(j) _ for some permutation vr. 
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It is worth noting that the standard Schwinger model is obtained as the special case 
A^P = Np = 1 and n+ = ?7-_, and the chiral Schwinger model corresponds to Np = Np = 1 
and n+ = 1 and n_ = (for review on previous work on these and similar models see 
Q). The class of models with only one photon field {Np = 1) and coupling constants enj^+ 
and erij- to the right- and left handed chiral components of the fermions, respectively 
(j = 1, . . . , A^F, and the nj^a are integers) was previously proposed and studied in |||]. For 
these models it is known that the gauge anomaly cancels whenever 

Np 

E(4.-4-) = o- (1) 

The simplest non-trivial (i.e. chirally asymmetric) solution is for A^f = 2, nj^+ = (3,4) 
and rij^- = (0,5), and therefore one sometimes refers to this model as the 3-4-5 model Q. 
Similarly one can find non-trivial examples for all Pythagorean Triple i.e. integer solutions 
of 0^ + 6^ = c^. 

For the general class of models with A^f fermions and A^p < A^f photons which we study 
we find the following conditions for gauge anomalies to be absent, 

Np 

E - nl.nl_) =0 Va, /3 = 1, . . . , ATp . (2) 

i=i 

A simple non-trivial example is for A'f = A^p = 2, nj = (3, 4), nj _ = (0, 5), = (—1,2), 
and rij _ = (2,1). To find solutions of these conditions for given A^p and as small A^f as 
possible seems to be an interesting generalization of the problem of finding Pythagorean 
Triples. 

The study of the 3-4-5 model in Ref. |jl| was in the path integral formalism. Here we 
use a Hamiltonian approach in the spirit of Ref. |p . This approach allows a rigorous con- 
struction of these model in the Hamiltonian framework using the quasi-free representation 
of boson- and fermion field algebras |^ and we also can solve the model using standard 
bosonization techniques |5|. Our general class of models (i.e. arbitrary A'f and A^p) is quite 
complicated, and it is somewhat surprising that it is possible to find the solution in such 
an explicit manner as we do in this paper. It is also intriguing to see the importance of 
the no- anomaly conditions in Eq. ^ at several different, seemingly unrelated points of our 
construction and solution. 

To simplify notation we explain our methods and computations in detail for the simplest 
case A^p = 1, and we are careful to do things such that the generalization to the case A^p > 1 
is easy. 

The plan of this paper is as follows. We first concentrate on the models with one photon 
field and an arbitrary number Np of fermion fields. After a formal definition we summarize 
the rigorous construction of these model in our framework. Anomalies are a consequence of 
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Schwinger terms which result from the normal ordering necessary to construct the fermion 
currents as well-defined operators on the Hilbert space of the model. Especially we compute 
the gauge anomalies from the commutators of the implementers of gauge transformations 
(= Gauss' law generators), and Eq. (||) is obtained as condition for a vanishing gauge 
anomaly. If and only if this latter condition holds a simple construction and solution of 
this model is obtained. As mentioned, these results for A^'p = 1 are presented such that the 
generalization to our general class of models with Np photons and N-p fermions is easy. The 
presentation of our results for this latter case are given in a rather short final paragraph. 

2. Notation 

Throughout this paper we consider Abelian gauge theories with massless fermions on two 
dimensional spacetime which is a cylinder, M x S"^, parametrized by coordinates {x'^) = 
{x^,x^) where x^ = t G R {= time) and x^ = x £ [—L/2, L/2] (= spatial coordinate; 
< L < oo is the spatial length). Our metric tensor is diag{l, —1). We use the Einstein 
summation convention only for spacetime indices iJ,,u = 0,1 but not for flavor- or spin 
indices. 

3. Formal definition of the model 

We now define in some detail the simplest non-trivial example for a chiral model without 
gauge anomalies. As described above, this model contains Np flavors of Dirac fermion fields 
'4>j,a, 'ipj,a and one Abelian gauge field here and in the following, a, a' = -|-, — are spin 
indices, and indices j,j' = 1, . . . ,Np distinguish the different fermion flavors. 
The model is formally defined by the Lagrangian 

C = + E ^J^^[-^9^. + e(n,,+P+ + n^- _P_)A^]V',- (3) 

i=i 

where F^^, = d^Ay - dyA^, and 

P± = i(l±73) (4) 

are chiral projections; here = (7^)^ct' are Dirac matrices which we take as 7*^ = cJi, 
7-*^ = — ia"2, and 73 = 7^7^ = CT3 (cTi.2,3 are the Pauli spin matrices as usual), and the 
real parameters Cj^a = enj^o- are coupling constants. The gauge group for this model is 
g = C°°(Rx 5^; U(l)) (= smooth U(l)-valued functions on spacetime), and the Lagrangian 
Eq. (^) is obviously invariant under the following gauge transformations, 

^Pj ^ (e^"^-.+^P+ + e^"^--^P_)V'j 

Af, ^ A^--d^x (5) 
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for all e'-*^ G Q. Note that the existence of the large gauge transformation e'-^^^'*^ = e'^/^ 
forces us to require that the integers (otherwise the large gauge transformations 

cannot be implemented in our model): The charges of the different fermion flavors have to 
be quantized. 

To motivate our construction of the model in the Hamiltonian framework below we 
recall some formulas from the formal canonical procedure; see e.g. Q. From the action Eq. 
(^) one computes the canonical momenta for the various fields of the model and obtains 
the following canonical (anti-) commutator relations[|'|^ 

[E{x),A,{y)] = -\6{x-y) (6) 
etc. where ip'^ = -07^ and E = Fqi. Moreover, with the notation 

Ho = dx'ip]{x)-fsi-idi)ijj{x) 

Ji^) = ^^nj,<T'4']{x)-i^Pailjj{x) 

N-e 

P{x) = rij^a^]{x)Pall)j{x) (7) 

the resulting Hamiltonian for the model can be written as 

H = Ho + J ^ dxi^^E{xf + eAi{x)J{x)j (8) 

and has to be supplemented by the constraint (= Gauss' law) 

G{x) = -diE{x) + ep{x) ~ . (9) 

We note in passing that we also have two Noether currents, a vector current (Jy) and an 
axial current (J^). These currents are given by 

Jy = /), Jy = J 

Jl = J, Ji = -p (10) 

and formally (i.e. prior to quantization) obey continuity equations, d^Jy = d^J^ = 0. We 
also introduce more general kinds of fermion currents 

Pj,a{x) = ipj^^{x)i;j^^{x) (11) 

which all are observables of interest for our model. 

^{o, b} = ab + ba and [a, b] = ab — ba 

^here and in the following we set t — and make explicit the dependence on the spatial coordinate only 
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4. Construction of the model 



We now outline a rigorous construction of this model using the representation theory of 
loop groups H. This construction amounts to representing the (Fourier modes) of the 
field operators wfa-, ^'^^ and the observable algebra of the model by closed operators 
on a Hilbert space TC such that the Hamiltonian H is represented by a self-adjoint operator 
on Tl. In this construction it is crucial to establish gauge- and Lorentz invariance. As we 
shall see, this will lead us to the condition in Eq. (|l|). 

The essential physical requirement determining the construction of the model and im- 
plying a non-trivial vacuum structure is positivity of the Hamiltonian H on the physical 
Hilbert space. As is well-known, this forces one to use a non-trivial representation of the field 
operators of the model. The essential simplification in (1+1) (and not possible in higher) 
dimensions is that one can use a quasi-free representation for the fermion field operators 
corresponding to "filling up the Dirac sea" associated with the free fermion Hamiltonian Hq, 
and for the photon operators one can use a naive boson representation. We now describe 
this representation in more detail. 

In the following it is convenient to work in Fourier space. We first introduce some useful 
notation which is such that in all equations the limit L — > oo is obvious. The Fourier space 
for even (periodic) and odd functions is 



A* = <i p = — n 



neZ} and A^^^ = <{ g = ^ ( n + ^ 



, , , n G Z 
L V 2/ 

27r 



respectively. For functions / on Fourier space we write Jj^*dpf{p) = J2peA* "r/(p) ^^'^ 

similarly for A*^^, and the corresponding 5-function satisfying / ^^dp' 6{p — p')f{p') = f{p) 
-La 

2tt "P,P' 



is 5{p — p') = -k-^p^p'- Our conventions for the Fourier transformed operators are 



^)[>{q) = / -= i(x)e - {q G A^,,) (12) 



(i.e. we use anti-periodic boundary conditions for the fermions), 

Mp) = pAiix)e-'P- (p G A*) (13) 

and in the other cases 

rL/2 

Y{p) = / dx Y{x)e-"P'' {p G A*) for Y = E,p, J, pj „ etc. (14) 

J-L/2 

Therefore the non-trivial canonical (anti-) commutator relations of the field operators be- 
come [Ai{p), E{p')] = i6{p + p') and 

{^i,f7(g),'0]/_ ^/(g')} = ^aySjj'S{q - q') 

{i^\%),i'\:l,iq')} = yp,p',j,j'. (15) 
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The model will be constructed in the full Hilbert space 7i = Wphoton "X" Wpermion- 
For 'Hphoton we take the Fock space generated by boson field operators (p) and (p) = 
^{pY 1 P £ (* is the Hilbert space adjoint) obeying the commutator relations 

[b{p), b\p')] = 5{p - p'), [b{p),b{p')] = yp,p' (16) 

and a normalized state Op G Wphoton such that 

b{p)np = Vp . (17) 

We then set 

ii(p) = i (bip) + 6t(-p)) E{p) = -| (bip) - b\-p)) (18) 

with a parameter s to be determined later. We recall that these requirements fix the Hilbert 
space "Hphoton completely. In this setting we now construct bilinears in the photon field 
using normal ordering : • • • : with respect to the state Q,p, e.g. '■b{p)b'^ {p') : = b^p')b{p) for 
all p,p' . 

For Hpcrmion we take the Fermion Fock space generated by operators il^j^fj^q) and 



'0jo-('?) = V'i,o-(9)* obeying the relations Eq. ([[SD and a normalized state Op E Wpermion 
such that 

V'+j(g)l^P = Vil,j(-g)l^P = V(7>0,Vi. (19) 

We note that the state Vt-p characterizing this representation of the fermion field algebra 
can be interpreted as Dirac sea associated with the free fermion Hamiltonian Hq. The 
presence of this Dirac sea requires normal-ordering : • • • : of the Fermion bilinears such 
as Hq = X^i /a* : 9 V'! (9)73^1 (^z) '■ and pj^a{p). This modifies the naive commutator 

^ odd ^ ' 

relations of these operators which follow from the canonical anti-commutator relations Eq. 
(p^) as Schwinger terms show up, see e.g. 0, ^. The relevant commutators for us are, 

[PjAP)^Pj',<T'{p')\ = <y&a,a'5jj'p5{p + p') 

[Ho,Pj,a{p)] = (^PPjAp) ■ (20) 

We also note that 

Pj,+{p)nF = Pj,-{-p)^F = Vp>0,Vi (21) 

which together with (|20|) shows that the pj^+{p) (resp. pj^^{p)) give a highest (resp. lowest) 
weight representation of the Heisenberg algebra. 

The (Fourier transformed) Gauss' law operators are now well-defined on 7i, 

G{p) = -ipE{p) + ep{p) (22) 
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with p{p) = J2j,a''^j,eTPj,a{p)- To determine a (possible) gauge anomaly we compute the 
commutators of these operators and obtain (we use Eq. (|20|)) 



[Gip),G{p')] = e^p6{p + p')J2 E • (23) 

j=l a=± 

We see that no gauge anomaly occurs if and only if the condition Eq. (||) holds. In the 
following we assume that this is the case. 

We now can give a precise meaning to the formal Hamiltonian H in Eq. (^) as follows, 

H = Ho + J^dp : (^^E{p)E{-p) + eAi{p)J{-p) + ttM^ Ai{p)Ai{-p)^ : . (24) 

Similarly as in the Schwinger model we have added a photon mass term (= last term 
on the r.h.s. of Eq. (p4[)) to restore gauge invariance which otherwise would be spoiled 
by Schwinger terms. Note that our normalization is such that this term formally equals 
iM2/f02dx^l(3;)^ i.e. M can be interpreted as photon mass. To determine the param- 
eter M'^ we compute {G{p),H\ and demand this to be always zero. By straightforward 
computation (we use Eq. (pO|)), 

2 

^^ = |^EE4.- (25) 

j = l(T=± 

We finally have to fix the parameter s in Eq. (|l^). We observe that H in Eq. ( p^ ) has the 
form H = iJo,F + -f^o.P + -f^int where Hq y = Hq and 

i^o,P = j ■■ (^^E{p)E{-p) + T,M^Ai{p)Ai{-p)^ : (26) 

are the free fermion- and Photon Hamiltonians. We fix s by requiring that -ffo,P has the 
form ] f^,dpu}{p)h\p)h{p). We thus obtain s = 2^/^T\M\ and uj{p) = |M|. Then i?o,F + -ffo,P 
obviously is defined as self-adjoint, positive operators on TC. Our results in this paper imply 
that also the interacting Hamiltonian H is well-defined on Ti.: the operator H is self-adjoint 
and bounded from below on TC. 

5. Bosonization and gauge fixing I 

As in the standard Schwinger model we can rewrite the Hamiltonians of our models as boson 
Hamiltonians by using the so-called Kronig identity Hq = |/^*dp Pj,'T{p)Pj,(7{—p) '■ 
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where : • • • : denotes normal ordering of the fermion currents (see for more details). We 
obtain 

Np 



H = J^dq : \l-E{p)E{-p) + E PjMPjA-p)j ■ 



(27) 



where we have defined 



Pj,aip) = Pj,a{p) + crenj^aAi{p) . (28) 

Note that [G{p), pj^a{p')] = i.e. the operators Pj,a{p) all are gauge invariant and thus 
observables for our model. 

The Hilbert space Ti. which we have obtained still contains gauge equivalent states i.e. 
states related to each other by static gauge transformations. To obtain the physical Hilbert 
space TCphys for our models we have to go through a procedure called fixing the gauge in 
the physics literature. We will do it in two steps. The first step will be done in this 
paragraph and eliminates the 'small' gauge transformations (i.e. the e'-^ G Q with x smooth 
and periodic functions on space 5^). This step is easy since the resulting Hilbert space 
^phys identified as the sub-Hilbert space of 7i which is annihilated by all the Gauss 

law operators. In a second step, which will be described in Paragraph 7 further below, we 
will account for large gauge transformations, too. 

For the first step we recall that the only gauge invariant degree of freedom of the Photon 

L /2 

field at fixed time is the holonomy /_£_y2 dx Ai{x). Due to the absence of a gauge anomaly 
we can therefore impose the gauge condition 

Ai{p)=5pflY (29) 

and solve the Gauss' law G{p) ^ (cf. Eq. (p^)) as 

E{p) = for p / 0. (30) 

ip 

This determines all components of E except the one conjugate to Y: 

^c' = ■ (^^' 

After that we are left with the {p = 0)-component of Gauss' law, viz. 

eQv - 0, Qv = m = Y.Y. ^i^PiA^) ■ (32) 

j=l a=± 

We thus can identify the Hilbert space of states invariant under all small gauge trans- 
formations as Wphys = L'^(R,dY) (g) Wpgj,jjjjojj where Hp^j.^-^^^ is the zero charge sector 
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of the fermionic Fock space. Moreover, by inserting Eqs. (p9[)-(|3l|), each gauge invari- 
ant operator on 7i becomes an operator on 'H'^Yiys- Especiahy the Hamiltonian becomes 
^ = X T.p>o ■ hp ■■ where 

^^ = -4^(^) + 2 ^ E (AAO) + -n.^^eYf (33) 

^ ^ j = l(T=± 

and 

■ hp-. = : \'^:;^p'^~P^p'yP^ ^1 ^ PjA-p)pjAp)j ■ ■ (34) 

This completes the first step of the gauge fixing procedure. As mentioned we have not yet 
accounted for the existence of the large gauge transformations. We will come back to this 
in Paragraph 7 below. 



6. Solution of the model 

The operators : hp : obviously all commute with each other. We thus can diagonalize the 
Hamiltonian H by diagonalizing each term : hp : separately. 

We first consider : hp : with p > 0. We introduce the boson operators 

-TnPj'Ap) forp>0 
Cj{p) = { V\P\ (35) 
■'^^ ' =Pj,~{p) for p < ^ ^ 



and c^j{p) = Cj{p)* = Cj{—p), so that [cj{p),c^j,{p')] = ^^j,j'Hp ~ p') Cj{p)^l = for all 
and p,p' ^ 0. We also find it convenient to define 

Cj{p) = Cj{p), CNp+j{p) = 4(p) i = 1, . . . A^F , (36) 

to fix p > 0, and suppress the argument p in the following. Then [cj,c^,] = ^Sjj'qj where 
qj = 1 and QNy+j = -1 for j = 1, . . . , Np. 

We then can write : hp : in matrix notation as : c^^ • he : where we defined = 
(c|, . . . , c\j^^) etc. and h is the 2Nf x 2Nf matrix 



e2 



+ n(g)n^ (37) 
2ttp 

where n-^ = (ni, . . . , n2Arp) with rij = nj^+ and n^^j^j = rij^^ for j = 1, . . . , Np (1 is the 
2Np X 2Np unit matrix here, and we use a standard tensor notation; the denotes the 
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matrix transpose). Defining hp = c"^ ■ he we get :hp: = hp— < ft, hpft > with 

(we used Eqs. (p5|) and @)). 

We now can diagonahze hp by finding a boson BogoUubov transformation C = Uc (U 
some 2A'^p x 2Np matrix) so that the operators Cj obey the same relations as the Cj and 
are such that hp = ■ DC = X]j=i ^j^l^j- It is easy to see that these two conditions are 
equivalent to 

UqU* = q and (U"^)*hU"^ = D (39) 

where q = diag{qi, . . . ,(?27Vf) and D = diag{u)i, . . . ,uJ2Np) (* and is the matrix adjun- 
gation and matrix inverse, respectively). To solve this somewhat unconventional diagonal- 
ization problem we note that Eq. (^) implies 

D^ = (qD)2 = U(qh)2U~^ , (40) 

and U is determined from this equation up to transformations U VU where V commutes 
with D^. We shah see that Eq. ( ^0[ ) corresponds to a standard diagonalization problem: 
the matrix (qh)^ is self-adjoint and thus can be diagonalized by a unitary matrix U. We 
thus can solve the problem in Eq. ( |39|) by first determining U and then making the ansatz 
U = VU with V commuting with D^. From Eq. ( [39| ) we then get the following condition, 

VUqhU*V"^ = qD (41) 



which again is a standard diagonalization problem and will allow us to determine V. 

We now compute D^ using Eq. (40). We write qh = pq+^eN^-^-i^eJ where we define 



ei = n/|n| (|n| = Vn^ • n) and CiVp+i := qei (we used |n| = |qn| and |np = M^). It is 
now crucial to observe that the condition Eq. (|l|) is equivalent to ■ qn = i.e. if there is 
no gauge anomaly the two vectors ei and e^^-^i are orthonormal. Moreover, in this case 
we can extend these vectors to a complete orthonormal real basis {ej}^^^ in R^^^ so that 
qej = e7Vp+j for j = 1, . . . , Ny- We thus obtain 

(qh)2 = pH + (ei ef + e^Vp+i ® 

from which we can immediately read off U and the matrix elements of D^ : denoting as Ej 
the standard basis in R^^f (i.e. (Ej)j/ = Sjj') we have 

U = ^ E, ej , (42) 



10 



and = + and for j = 1, A^f + 1 and = otherwise. We then compute 



M. 



UqhU* = ^Etvp+j ® eJ + pEN^+j eJ + pEj ® E^p+^. (43) 



P 



with Ml = M and Mj^i = 0, which shows that we can determine V by diagonahzing 2x2 

/ p\ 
matrices of the form a/2 . We find 



2Afp 

V = ^ Ej FJ (44) 



where 
Thus 

We thus obtain 



pE,-±>2 + M/E;Vp+, ^ ^ 

'i,Wp+i = , ^ „ . • (45) 



00 j = -WTVp+i = V^J' + for j = 1, . . . , iVp . (46) 



TVp 

^ = E 7^^+^ (C^.(P)^C^.(P) + C,(-p)C,(-p)t) (47) 
i=i 

where we used CN^+jip) = C'j(— p). This completes the diagonahzation of the operators 
-.hp-. . 

We are left to diagonalize Hq. We note that this is the part of the Hamiltonian which 
contains the operators Y and pj^ai^) which are not invariant under the large gauge trans- 
formation e'^/^ but transform as follows, 

Y ^ V-i 

e 

PjA^) /Oj,a(0) . (48) 

It is therefore not immediately obvious that ho is indeed invariant under the large gauge 
transformations. To make this invariance manifest we note that the operator Qa = 
^j^„(ynj_fjf)j^(j(fi) changes under the transformations Eq. (4S) as Q\ Qa + |np, and 
therefore 

Y :=¥+ -t^Qa (49) 

e nr 
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is indeed invariant. By straightforward computations we obtain 

\M\ 



^0 = ^ (Ci(O)tCi(O) + Ci(0)Ci(0)t) +C (50) 



where ^{^^(0) =V ^(^|M|)-i/2^ + (^|M|)V2 y and 



\{j:i:hM'-T^,Qi\ (51) 




2 \j = la=± 

ah are invariant under the transformations Eq. (^8|). We also have [Ci(0), Cj(0)] = ^ and 
thus see that /iq is essentially a harmonic oscillator Hamiltonian. 
We can combine our results above in a compact form as follows, 

^ = E y ^}p^p' + M]c]{p)C,{p) + —C + £o, (52) 

with the constant 

^o = E(^ + :^/ dp{Jp' + Mf-{p+:^)]] (53) 




(note that actually only the term with j = 1 contributes to i?o)- This shows explicitely that 
the Hamiltonian H is self-adjoint and bounded from below, and since C is non-negative 
(we will show this further below) the constant Sq is equal to the the ground state energy. 
Moreover, we also see explicitly that our model has a relativistic spectrum: the physical 
degrees of freedom correspond to one massive- and Np — 1 massless boson fields. 

We now briefly describe how to construct a groundstate for our model. We note that 
with the explicit formulas given above it is straightforward to construct a unitary operator 
U{p) in terms of the operators Cj\p) such that Cj{p) = U{p)cj{p)U{p)* . Moreover, one 
can check that U = Y{p^oU{p) defines a unitary operator on Wpgj.jj^joj^. It is then easy to 
see that the state lpq{Y)UQ.y G ^phys '^^^'^ '/'o(^) oc exp(— 7r|M|^y^) is annihilated by all 
operators Cj{p) and C and thus a ground state of H. We thus have found one groundstate 
for our model. However, this state is highly degenerate, and it is actually not invariant 
under the large gauge transformations Eq. (|48|). In the next paragraph we will discuss this 
issue in more detail. 

We finally mention that in a complete solution of the model one also needs to find 
the Green's functions i.e. vacuum expectation values of gauge invariant combinations of the 
fermion- and photon field operators. All these Green's functions can be computed explicitly 
by using the so-called boson-fermion correspondence which allows to write the fermion field 
operators in terms of the fermion currents. These computations are similar to the ones for 
the usual Schwinger model (see e.g. Q) but beyond the scope of the present paper. 
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7. Gauge fixing II: Vacuum structure and all that 

We now describe the structure of the Hilbert space of our model and then perform the 
second step of the gauge fixing procedure described already above. 

We first recall the well-known structure of the fermion Fock space TYFermion^ this space 
can be generated by the fermion current pj^a-{p), P 7^ 0, together with unitary operators 
Rj^fT which obey the relations 

R'^Pji^a' {p)Rj,a = Pj',a' {p) + 5pfi5jjia5„y (54) 

and which interpolate between different sectors labeled by the eigenvalues of the charge 
operators. Thus for all m = (?ni^+, . . . , m7Vp,+ , . . . ,m7Vp,-) G 7?^^ , the operators 

= R^x^ . . . rZ!+^rt:-- ■ ■ ■ RZr (55) 

commute with all /ip>o and U, and if ^' is a vector in ^Fermion with pj^o-(0)\I' = Vj, a then 
Pj^a{0)IV^^ = rn-j^o-R-™^) which implies that R™^' is an eigenvector of Qy and Qa with 
eigenvalues n-^ • m and n-^ • qm, respectively. This implies that if such a vector ^' is also 
an eigenstate of all /ip>o and if € -^^^(R) an eigenstate of /iq then all states 

V9(y)R™^ = ip(Y+ "^^^ ^ R'"* 

with n-^ • m = (= charge zero condition) are eigenstates of H with eigenvalues which 
are of the form £ + c(m) and differ only by the contribution from C Eq. (|5l|), c(m) = 
^^^■2 (|mp|np — (n-^ -qm)^). Note that the latter is always non-negative due to the Cauchy- 
Schwartz inequality. We thus see that all these states with m = fcqn {k integer) are 
degenerate, and especially all states ipo (y + |^ H'^^'^UQ are groundstates for our models. 
This degeneracy actually is explained by the invariance under large gauge transformation 
Eq. (^): This transformation acts on states in Ti'^^y, as ip{Y)^ ^ 99 (y + i) Ri"^' . 

We now come to the second step of our gauge fixing procedure. Our discussion above 
implies that the states 

iffiY) = e'^V (y + -) R^''^"* (56) 



(9 real) have simple transformation properties under large gauge transformation Eq. (^), 
ifl^ — > e^^[^, ip]^ . We thus can define Wphys as the vector space spanned by all [^', p]^{Y). 
However, this does not yet make Wphys into a Hilbert space: a simple computation shows 
that 

< [^l,ipi]^, [*2, f2f >Wphys= - G') < V'l], [*2, ^2] >phys 
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where < [^i,ipi],[^2,^2] >phys=< *i,^2 >WF„niion< ^i^^2 >l^(r) is independent of 6 
and 6'. We see that the vectors [^,ip]^{Y) are actually not contained in Wpj^yg (they do 
not have finite norm). The remedy of this problem is a simple multiplicative regularization 
i.e. 'dropping the infinite constant 5(0)'. This is equivalent to using < •,• >phys as inner 
product in Wphys which is well-defined. This completes the construction of the model. 

8. Generalization to an arbitrary number of photons 

We now describe how the above results generalize to a large class of models with A^p different 
photon fields A'^ where a = 1, . . . , A'p. These models are (formally) given by , 

I Afp Afp Afp 

^ = "4 E F^uiFn^" + E ^J^'i-^d, + ^ e^{nl^P+ + n-_P„)A-]V^, (57) 

a=l j=l a=l 

where = f?^A" — d^A'^ and the charge units corresponding to the different gauge 
fields can be different. This model obviously is invariant under transformations belonging 
to the gauge group G = C°°(R x 5^;U(1)^p), A'^ ^ A'^ - etc., and as before the 

existence of large gauge transformations requires all the n"^ to be integers. The canonical 
procedure and the construction of the models with A'^p = 1 generalize with minor changes: 
now we have A^p copies of the photon fields and correspondingly A^p copies of the Gauss' 
law operators etc. To be more specific: We now have [Af (p) , (p')] = iS'^'^S{p + p'), a 
Hamiltonian 

H = Ho+ dpJ2--\ -rE%p)E'^{-p) + e"AUp)r{-p) + vrif (p) T"/^i^^(-p) : , 

(58) 

and Gauss' law operators G"'{p) = —ipE°'{p) + p"'{p) where 

P"(P) = E <<rPjAp) > J'^iP) = E ^^JUp) (59) 

with pj^u{p) as before. The model has no gauge anomalies if and only if all commutators 
[G"(p), G'^(p')] vanish, and similarly as for A'^p = 1 we obtain the conditions in Eq. (^). Note 
that the mass term we have to add to the naive Hamiltonian depends on a A^p x A'^p-matrix 
T"'^ which is determined such that H commutes with all G"{p). We thus obtain 

a /3 Np 

^"' = ^EE-"X.- (60) 

j=l a=± 
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We note already here that T = (T"'^)^^^-^ is a self-adjoint, real, non-negative A^p x A'^p 
matrix, and therefore we can write 

T°^ = ^M2(a^)"(a^)/3 (61) 

7=1 

where (a^)° are the components of the orthonormal eigenvectors of T and the corre- 
sponding eigenvalues. For later convenience we assume that rank(T) = A'p i.e. that all the 
are non-zero.0 

The generalization of the representation Eq. ( p!8| ) of the photon fields etc. generalizes 
in a straightforward manner (we only mention that the generalization of the free Photon 
Hamiltonian now becomes Hq^p = J2a I A*'^P\^^a\b}^{p)ba{p) where are the eigenvalues 
of the matrix T). Moreover, also the generalization of Eq. ( p7[ ) is obvious where the gauge 
invariant currents now are Pj,a{p) = Pj,a{p) + ^J2a=i^°'''^'j,a'^iip) (note that due to Eq. 
(^ these latter currents indeed commute with all Gauss law operators). 

We now come to the solution of the models: the gauge fixing condition generaliz- 
ing the one in Eq. (^) obviously is Ai{p) = 5pfiY°'. Imposing that condition, the 
Hilbert space of states invariant under all small gauge transformations becomes Wpj^yg — 
-L^(R^p, dY^ ■ ■ ■ dy^P) 8 Wpgj.jjjjQj^ where the zero charge sector in the fermion Fock space 
now is defined such that 

= E E ^.Vj--(O) = Va = l,...,iVp (62) 

j=l a=± 

on Wpgrmion (again, these latter conditions come from the {p = 0)-components of Gauss' 
law). Moreover, the Hamiltonian can again be written as H = ^J2p>o '-^p- ■ 

For p > we can use the matrix notation introduced in Paragraph 6 and write hp = 
• he where h = pi + with 

JVp Q, (J 

Y = y ® (n")^ (63) 

27r 

and (n°)"^ = {nf _^_, . . . , . . . , _) etc. as before. We now show how to diag- 

onalize these Hamiltonians, following the method explained in Paragraph 6: we compute 
(qh)^ = p^l + Y + qYq where we used YqY = which follows from Eq. (^). We now 
observe that Y is a self-adjoint, real, non-negative matrix, and it therefore can be written 
as 

Y = E ^j^i ® 4 (64) 



'We believe that this assumption could be easily dropped. 
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where the e,,- are orthonormal and > Ml > • • • > • "^^^ ^iid ej can computed 
by diagonalizing the A^f x -^f matrix Y. We also observe that the non-zero eigenvalues Mj 
of the matrix Y are identical with the non-zero eigenvalues of the matrix T = (T"^)^^^^ 
defined in Eq. (|60|). Thus only the Mj with j = 1, . . . , rank(T) = A'p are non-zero. (To see 
this, note that if (flj)" are the components of an eigenvector of T with non-zero eigenvalue 
MJ, then fj = X^a is an eigenvector of the matrix Y with the same eigenvalue 

Mj. One can also check easily that the ij span a space of dimension equal to the rank of 
the matrix T which we assumed to be equal to A'^p.) 

Defining e^^j^j = qej we obtain a complete orthonormal basis in R^^^ (orthogonality 
of the Gjxtvf ^j>Np again follows from Eq. (|2|)), and we can write 



Afp 



(qh)2 = / l + g m| (ej ej + ejVp 



Thus U Eq. ( ^2|) diagonalizes the matrix (qh)^, and the eigenvalues of this latter matrix 
are oj] = p'^ + Mf where Mj^^^j = Mj. It is then easy to check that all the Eqs. (|3|)-(|4|) 
remain true also in the present case, and one finally obtains a representation of hp as in 
Eq. (P). 

We now turn to /iq- After some computations we obtain, 
where Y^ = Y^ + Ef^i^'^r^Qi, Qa = Ej,a<^rilJj^^{0), and 

c = M E E p^Aof - E qu-^-Wa ■ (66) 

Note that Y", C and /iq all are invariant under the large gauge transformations 

Y'^ > Y'^ — -— 

Np 

PjAO) ^ /5,,.(0) + E K,-^" (67) 

a=l 

where k= {k^,...,k^p) G Z^p. 

Introducing Zj = J2a{0'j)"^"' where (a^)" the components of the eigenvectors of the 

matrix T (cf. Eq. ^)) and 0^(0) =+ ^(vr|Mj |)-V2_|_ + (^|jv/^.|)i/2 we can write 

Afp I ^ I 

^0 = E V {cm^C,{0) + C,(0)C,(0)t) + C . (68) 
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Again we can combine our results and write H as in Eqs. (p2|)~(|5^), and we see explicitly 
that our model has a relativistic spectrum: we have A^p massive and {Np — Np) massless 
bosons. 

It is straightforward to extend the construction of the unitary operator U on Ti-'p^^^^^^ 
diagonalizing all the hp and then check that ipQ{Y^, . . . ,Y^P)L{Qp G ^phys with ipo = 
exp(-^ Ea,/3,7 \M^\ia'y)''ia'yfY°'Y'^) (cf. Eq. (|l|)) is one groundstate of the model. One 
then can check that for any ip G L^(R^p) which is an eigenstate of /iq and any ip G '^Fermion 
which is a common eigenvector of all hp, the state 



V{Y\...,Y 



Nr- 



with (n°)"^ • m = 1 for all a, is an eigenstate of H with an eigenvalue of the form E + c(m) 
where 

c(m) = — - ^(ej • qm)^ 
^ ^ i=i 



with ej the orthonormal eigenvectors of the matrix Y in Eq. (^) . With that one can check 
that the eigenstate of H which also are invariant under the large gauge transformations are 

^^""^ kez^p 

X,^ l^yl + ^, . . . , y^P + . . . Rfc-Pqn^P ^ (gQ) 

where now we have A^p real 9 parameters. Similarly as for A'^p = 1 the physical Hilbert space 
Wphys of the model is spanned by the states (f]^^'---'^\Y), and one needs to renormalize 
the inner product of these states, i.e. 'drop the infinite constant d{0)^^\ to get a proper 
inner product on Wphys- 
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